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Motivated by the recent experiment at ENS [ V. Bretin, S. Stock, Y. Seurin and, J. Dalibard, 
Phys. Rev. Lett. 92, 050403 (2004)], we study a rotating (non-)interacting atomic Bose-Einstein 
condensate confined in a harmonic-plus-gaussian laser trap potential. By adjusting the amplitude 
of the gaussian laser potential, one can make quadratic-plus-quartic potential, purely quartic poten- 
tial, and quartic-minus-quadratic potential. We show that an interacting Bose-Einstein condensate 
confined in a harmonic-plus-gaussian laser trap breaks the rotational symmetry of the Hamiltonian 
when rotational frequency is greater than one-half of the lowest energy surface mode frequency. We 
also show that by increasing the amplitude of the gaussian laser trap, a vortex appears in a slowly 
rotating Bose-Einstein condensate. Moreover, one can also create a vortex in a slowly rotating 
non-interacting Bose-Einstein condensate confined in harmonic-plus-gaussian laser potential. 
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The rotation of a macroscopic quantum fluid exhibits 
interesting counter-intuitive phenomena. For example, 
an atomic interacting(repulsive) Bose-Einstein conden- 
sate (BEG) confined in a rotating harmonic trap pro- 
duces quantized vortices for a sufficiently large rotation 
[1-3]. The theoretical calculation for flc based on purely 
thermodynamic arguments is significantly smaller than 
the observed value of which is the order of ~ O.Twq, 
where ujq is the radial trap frequency. The vortex for- 
mation of a harmonically trapped atomic BEG is related 
to the dynamical instability of the lowest energy surface 
mode excitations whose energy scale is set by the har- 
monic potential [4]. Above a certain rotation frequency 
of a harmonically symmetric trapped interacting BEG, 
the system itself starts deforming from a circular shape 
to an elliptic shape and hence it breaks the original ro- 
tational symmetry of the Hamiltonian. This has been 
predicted only for a harmonically trapped BEG [5] and 
detected at the ENS experiment which leads to the vor- 
tex nucleation [6] . The vortex nucleation starts when the 
average angular momentum of each particle is one which 
has been measured experimentally [2]. 

The harmonically trapped BEG becomes singular 
when the rotation frequency is equal or greater than the 
harmonic trap frequency, because the outward centrifu- 
gal force counteracts the inward force from the harmonic 
trap. One can eliminate the singularity at 51 > 1 by con- 
sidering an additional stiffer radial potential (say, quar- 
tic potential for simplicity). There is a growing inter- 
est about the effect of an anharmonic potential on the 
properties of a rotating BEG [7-13]. In spite of the fact 
that the fl^ in a harmonically trapped BEG calculated 
from the thermodynamic arguments does not match with 
the experimental values, many authors have studied Qc 
in a BEG confined in quadratic-plus-quartic trap based 
on the thermodynamic arguments [8-11]. In this work, 



we show that a vortex appears in a BEG confined in a 
quadratic-plus-gaussian laser potential due to the spon- 
taneous shape deformation of the system as it is happens 
in a BEG confined by harmonic trap only and calculate 
the correct critical rotational frequency to create a single 
vortex. We also show that by increasing the magnitude of 
the laser trap, a vortex can nucleate in a slowly rotating 
Bose gas. 

Recently, the quadratic-plus-quartic potential has been 
achieved experimentally by superimposing a blue de- 
tuned laser beam to the magnetic trap holding the atoms 
[15]. The effective external potential is 
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where e = 



f is the deformation parameter due to the 



stirring potential which rotates the system and creates an 
anisotropic potential in the xy plane. Also, 2ujq — uj'^+ujy 
and the potential created by the gaussian laser beam is 



U{r) = Uoe- 
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At the ENS experimental set up [15], the laser waist 
is w — 25/im and the amplitude of the laser beam is 
Uo ^ 1.242 X 10""^° J. The harmonic trap frequencies 
are = 27r x 75.5Hz and ojz — 2% x llHz. Typically, 
the size of a condensate is the order of few iim. Since 
< 1, one can expand the gaussian laser potential. 
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The resulting potential can be written as. 
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where the co-ordinates f, i, y, and z are in units of the 
harmonic oscillator length ao = = 1.475/im. 

Here, k = and A = k(^)^. For the given val- 

ues of the w, Ua and = 27r x 75.5Hz, k = 0.24 and 
A = 8.297 X 10"" at the ENS experiment [15]. The value 
of k renormalize the oscillator frequency {loq) of the mag- 
netic trap and also the anharmonic term A depends on 
the value of k. At the ENS experiment, k is small. In 
this paper, we are interested to study when k is large. 
When fc < 1, fc = 1 and A; > 1, the effective potential is 
quadratic-plus-quartic potential, only quartic potential 
and quartic-minus-quadratic potential, respectively. If 
k > 1, the effective external potential becomes Mexican 
hat structure. In principle, fc > 1 can be easily obtain 
in the current experimental set up at ENS by increasing 
Uo from Uq = 1.242 x 10"^° J to Uq > 5 x 10-^° J. For 
simplicity, we study quasi-2D system since the elongated 
condensate expands radially and contracts axially when 
it rotates rapidly. The equation of motion of the conden- 
sate wave function 'ijj{r) is described by the mean- field 
Gross-Pitaevskii equation. 
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Here, 52 = 2\/2-K%LOza'a'z [16] is the strength of the mean- 
field interaction, a is the s-wave scattering length and 

= Also, Lz = xpy — ypx is the ^-component 

of the angular momentum operator and CIq is the trap 
rotation frequency. When e = 0, it preserves the circular 
symmetry of the system. 

One can write down the Lagrangian density corre- 
sponding to the quasi-2D system as follows: 



Here, we use the time-dependent variational method 
to study the properties of a rotating BEG confined in 
a quadratic-plus-laser potential. In order to obtain the 
evolution of the condensate we assume the most general 
Gaussian wave function, 



ij{x,y,t) = C{t)e-^ 



^\a{t)x^+/3(t)y''-2i(t)xy] 



(5) 



where C{t) is the normalization constant. Further, 
a{t) = ai{t) + ia2{t), I3{t) = Pi{t) + iP2{t) and -f{t) = 
71 {t) + ij2 {t) are the time-dependent dimensionless com- 
plex variational parameters. The ai and /3i are inverse 
square of the condensate widths in x and y directions, 
respectively. The above mentioned order parameter de- 
scribes only the vortex free condensate (irrotational sys- 
tem) as the phase {S{x,y) = jxy) corresponds to the 
irrotational velocity flow. 



We obtain the variational Lagrangian L by substitut- 
ing Eq. (5) into Eq. (4) and integrating the Lagrangian 
density over the space co-ordinates, 

^ = l^[-(/3ia2 + ai/32- 27172) + (Qi+/3i)I> 



- 2(a2 + /32)7i72 + (1 - k){ai + /3i) + e(/3i - ai) 
A (3aj + 2ai/?i + 3/?f + 47^) 
2 D 
+ PD^/^ + 20(71 (/32 - a2) + 72(^1 - A))], (6) 



+ 



where D = \J a\(5\ — ^\ and P = 2y f ''^^ 
mensionless parameter that indicates the strength of the 
mean-field interaction and it can be positive or negative 
depending on the sign of the s-wave scattering length o 
[14]. 

The variational energy of the rotating condensate at 
equilibrium is given in terms of the inverse square width 
of the condensate along the x and y directions and the 
phase parameter ^72, 
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where Q = ilo/u^o- We have used the fact that a2 = (^2 = 
7i = to obtain the above equation. 

One can get the equilibrium value of the variational pa- 
rameters, aio = X, /3io = Y and 720 = by minimizing 
the energy with respect to the variational parameters. 
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Eqs. (8), (9) and (10) describes how the shape of the 
condensate changes due to the change of the amplitude 
of the laser beam as well as external rotation. 

The average angular momentum per particle is given 
by, 
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where 77 = is the ratio of the square of the widths 
along the x and y directions. Eq. (11) exphcitly shows 
how the angular momentum is transfered to the trapped 
BEC due to the shape deformation with the rotation fre- 
quency. The vortex would appear when < L^, >= Nh. 
Therefore, one can estimate the critical angular frequency 
to create a single vortex from the relation: < Lz >= NU. 

First we calculate the quachiipolc mode frequency of 
BEC confined in a magnetic-plus-laser potential. We ex- 
pand the Lagrangian in the following way: a = X + 6ai, 
P = Y + 6 Pi, and 7 = 571. We keep only the second 
order deviations from their equilibrium values. Then we 
calculate the Lagrangian quadratic in the deviations and 
using the Euler-Lagrangian equation of motion, we get 
the following two coupled equations of 6ai and 6(3i: 
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For an isotropic trap, e = 0, and we set 6ai = —6/3i to 
calculate the lowest energy quadrupole mode frequency 
which is given by, 




FIG. 1. Plots of quadrupole mode frequeney coq vs inter- 
action parameter P for fc = 0.9, A = 0.0031(upper curve), 
k = 1.0, A = 0.0034(middle curve) and k = 1.1, A = 0.0038 
(lower curve). 

Fig.l shows that when fc < 1, the ujq decreases as P 
increases and after a critical value of P, uiq start increas- 
ing slowly. On the other hands, for /c > 1, increases 
as P increases. 
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FIG. 2. Plots of quadrupole mode frequency uiq vs k for 
P = O(dot-dashed), 10 (dashed), 200 (solid). 
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where Rq is the equilibrium radius of the system, and 
this can be obtained from the real solution of the cu- 
bic equation: (1 + f )i?^ - (1 - k)Ro - 4A = 0. The 
above quadrupole frequency ujq is valid for all interaction = —2 is 

(repulsive, attractive) strength. The quadrupole mode 
frequency ujq vs P for different values of fc, i.e. for dif- 
ferent configuration of the efi^ective potential, is shown in 
Fig.l. 



The quadrupole mode frequency iVq vs k for different 
(14) values of P is shown in Fig. 2. For fc = = P, it re- 
produces the known result for the quadrupole mode fre- 
quency of a harmonically trapped noninteracting BEC. 
With the increasing of fc, Uq of an ideal Bose gases is 
decreasing faster than the interacting cases. 

The centrifugal term —ClLz shifts the quadrupole mode 
frequency by —20. The quadrupole mode frequency with 



iV-2 = W 4(1 - fc) + ^ - PRl - 2fl. (15) 
V ^0 

We will discuss the effect of the harmonic- plus-gaussian 
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laser trap potential on the possibility of a vortex forma- 
tion in a slowly rotating BEC. Using the three coupled 
equations (8), (9) and (10), we get the following two cou- 
pled polynomial equations: 
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(fi^ - 1 + fc - e) - X{bX^Y + 7XY^ + X'^ + W^) (16) 
and 

= (X3y3 + 2x2y4 + xY^){i + ^s^^) 

- XF^(30^ + 1 - A; + e) + {X^Y + 2X^y^) x 

(f2^ - 1 + fc - e) - A(5Xy2 + IX'^Y + iX^ + Y^). (17) 



the system is destabilized. Therefore, there is no spon- 
taneous shape deformation and hence a vortex can not 
appear in the noninteracting gas confined by harmonic 
trap only. When k ^ Q, the the lowest energy surface 
mode frequency becomes less than 2wo and an instability 
will occur when O < wq at which the system is stabilized. 
There is a spontaneous shape deformation and one can 
produce a vortex in a slowly rotating noninteracting Bose 
gases confined in a harmonic-plus-gaussian laser trap. 



For simplicity, we put e = and solve those two coupled 

equations numerically. We find that spontaneous defor- 
mation occurs when the rotational frequency is greater 
than one-half of the quadrupole mode frequency Wg, and 
start transferring the angular momentum to the system 
[see Eq. (11)]. This circular symmetry breaking is due to 
the tendency against an instability of the lowest energy 
surface mode frequency with = —2. The variation of 
rj with the rotational frequency f2 for e = and k = 0.9 
is shown in Fig. 3. The same analysis can also be done 
for other values of k and P. 
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FIG. 3. The variation of rj vs the rotational frequency 
for P = 200, k = 0.9, and e = 0. 

For k = (only harmonic potential) and P is large, 
one can produce a vortex when f2 ~ 0.7. For fc ^ 0, 
the critical rotational frequency (flc) to create a single 
vortex is small (fie < 0.7) compared to harmonic trap 
case {Qc ~ 0.7). By adjusting the strength of the gaus- 
sian laser beam, one can produce a vortex in a slowly 
rotating BEC. For a purely harmonic trapped ideal Bose 
gases, the lowest energy surface mode frequency is 2a;o 
and hence instability will occur when O = at which 



FIG. 4. Critical rotational frequency fie 
' = 200(sohd line) and P = (dashed line). 



vs k for 



Fig. 4 shows how fic is decreasing with the parameter 
k for different values of P. As k increases, Oc is also de- 
creasing. Hence there is a possibility of vortex formation 
of a slowly rotating BEC if k is finite. When k < 0.8, 
flc is small in the Thomas-Fermi regime compared to the 
non-interacting case. When k > 0.8, Q,, is large in the 
Thomas-Fermi regime compared to the non-interacting 
case. It implies that the two-body interaction helps to 
create a vortex when k < 0.8, but the two-body interac- 
tion do not help when k > 0.8. 

In this work, we have shown that by adjusting 
amplitude of the gaussian laser beam one can make 
quadratic-plus-quartic potential, purely quartic poten- 
tial and quartic-minus-quadratic potential. We have also 
shown that an interacting BEC confined in a harmonic- 
plus-gaussian laser tarp breaks the rotational symmetry 
of the Hamiltonian due to the instability of the lowest 
energy surface mode. It is argued that by increasing the 
amplitude of the gaussian laser trap, a vortex can be 
created in a slowly rotating interacting as well as nonin- 
teracting BECs. 
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